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We report the realization of quantum microwave circuits using hybrid superconductor- 
semiconductor Josephson elements comprised of In As nanowires contacted by NbTiN. Capacitively- 
shunted single elements behave as transmon qubits with electrically tunable transition frequencies. 
Two-element circuits also exhibit transmon-like behavior near zero applied flux, but behave as flux 
qubits at half the flux quantum, where non-sinusoidal current-phase relations in the elements produce 
a double-well Josephson potential. These hybrid Josephson elements are promising for applications 
requiring microwave superconducting circuits operating in magnetic field. 

PACS numbers: 74.45.+c, 62.23.Hj, 84.40.Dc, 85.25.Hv 


In superconducting circuits, macroscopic degrees of 
freedom like currents and voltages can exhibit quantum 
mechanical behavior. These circuits can behave as arti¬ 
ficial atoms with discrete, anharmonic levels whose tran¬ 
sitions can be driven coherently [1]. In the field of cir¬ 
cuit quantum electrodynamics (cQED), these artificial 
atoms are coupled to resonators to perform microwave 
quantum optics in the solid state [2, 3]. Over the last 
decade, cQED has also grown into a promising platform 
for quantum information processing, wherein the ground 
and first-excited levels of each atom serve as an effec¬ 
tive qubit [4]. To date, implementations of supercon¬ 
ducting quantum circuits have relied almost exclusively 
on aluminum/aluminum-oxide/aluminum (A1 /AlOx/Al) 
tunnel junctions as the source of non-linearity without 
dissipation. However, many exciting applications require 
magnetic fields 0.5 T) at which superconductivity 
in aluminum is destroyed, calling for an alternative ap¬ 
proach to realizing microwave artificial atoms. 

Recent advances in materials development and 
nanowire (NW) growth have enabled the development 
of superconductor-semiconductor (super-semi) structures 
supporting coherent charge transport without dissipa¬ 
tion [5], and providing signatures of Majorana bound 
states [6]. To date, super-semi-super Josephson elements 
(JEs) have been studied exclusively in quasi-DC trans¬ 
port [7-10]. Building microwave circuits operating in the 
quantum regime, in which transition energies between 
levels exceed the thermal energy, offers news ways to in¬ 
vestigate the physics of hybrid super-semi structures us¬ 
ing spectroscopy [11-13]. 

In this Letter, we report the realization of microwave- 
frequency cQED circuits made from hybrid JEs based on 
In As NWs contacted by NbTiN. Capacitively shunted 
single JEs behave as weakly anharmonic oscillators, or 
transmons [14], with transition frequencies tuneable by 
the field effect, i.e., voltage on a proximal side gate. 


Double-element devices show similar transmon-like be¬ 
havior at zero applied flux, but behave as flux qubits [15] 
near full frustration owing to a double-well Josephson 
potential arising from non-sinusoidal current-phase re¬ 
lations (CPRs). We observe microwave driven transi¬ 
tions between states with oppositely flowing persistent 
currents, manifesting macroscopic quantum coherence. 
Their fabrication from magnetic-field compatible mate¬ 
rials makes these JEs promising for applications requir¬ 
ing quantum microwave circuits withstanding magnetic 
fields up to 0.5 T, such as the braiding of Majorana parti¬ 
cles and the interfacing with coherent quantum memories 
based on polarized electron-spin ensembles. 

Our chip [Fig. 1(a)] contains multiple capacitively 
shunted single and double NW JEs coupled to dedicated 
transmission-line resonators for control and readout us¬ 
ing a common feedline [Fig. 1(b)]. The chip contains side 
gates for electrostatic tuning of some single-junction de¬ 
vices and current-bias lines for threading flux through 
the loops of split-junction devices. We created each 
JE by deterministically placing an InAs NW between 
the leads of a pre-patterned NbTiN interdigitated ca¬ 
pacitor (IDC) and contacting the NW to each lead in 
a subsequent NbTiN deposition. The charging energy 
Eq ~ h x 300 MHz of the devices is chosen much smaller 
than the estimated Josephson coupling energy Ej of the 
NW junction, as in conventional transmon devices [14], 
leading to a weakly anharmonic energy spectrum (en¬ 
ergies Ei) of circuit plasma modes. We first verify the 
presence of the non-linear NW circuit by measuring the 
feedline transmission near the fundamental frequency of 
the coupled resonator [Fig. 1(d)]. The Jaynes-Cummings 
interaction leads to different resonator frequencies f c and 
/bare at single- and many-photon probe levels, respec¬ 
tively [16, 17]. We then search for the qubit transition 
frequency /oi of the NW circuit by monitoring feedline 
transmission at f c while sweeping a second tone [18] near 
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FIG. 1: Hybrid Josephson elements in cQED. (a) Overview of 
cQED chip allowing control and readout of NW circuits using 
dedicated resonators (green) coupled to a common feedline 
(blue). For readout, a microwave tone with frequency f c is 
applied near the fundamental of the resonator coupling to the 
NW circuit under study. The signal is amplified and downcon- 
verted to 1 MHz using a local oscillator at /lo for subsequent 
digitization and processing. Additional controls on a subset 
of devices include side gates (orange) for electrostatic tuning 
of carrier density in the NW of single-junction devices, and 
short-circuited transmission lines (purple) for threading flux 
through the loops of split-junction devices. Measurements on 
devices 1-3 are described in the main text. See table SI for an 
overview of the other devices, (b) Optical zoom-in of Device 
1, containing a single-junction NW circuit (red), (c) Scan¬ 
ning electron microscope (SEM) image of an InAs NW (blue) 
contacted by NbTiN electrodes (red) separated by 500 nm. 
(d) Normalized feedline transmission as a function of readout 
power. The resonator shifts from f c — 3.9464 GHz at single¬ 
photon level to /bare = 3.9470 GHz above 10 5 photons. This 
shift confirms the coupling of the resonator to a non-linear 
circuit. 


the estimated frequency /bare+(/bare - fc) / 9 2 , where g is 
the coupling strength between NW circuit and resonator. 

We first investigate the electric-field effect on the NW- 
circuit spectrum [Fig. 2(a)]. Device 2 has one NW 
junction (measured length L ~ 550 nm) and a proxi¬ 
mal side-gate electrode for tuning the carrier density in 
the NW. We observe fully reproducible fluctuations [19] 
in /oi as a function of the side-gate voltage V g , indi¬ 
cating phase-coherent diffusive charge transport in the 
NW [5, 20]. Using the plasma-oscillation relation [14], 
Ej ~ f^/SEc, we determine the root-mean-square 
Josephson-energy fluctuation yj (S 2 Ej)/h ~ 2 GHz in the 
V g = 0 — 10 V range. Matching this scale to the Thouless 
energy [21], Et h = hD/L 2 , and assuming highly trans¬ 
parent contacts [9], we estimate the diffusion constant 
D « 40 cm 2 /s. This value is typical for InAs wires [7]. 

Side-gate tuning of the NW junction offers a new 
means to control the spectrum of transmons. Decreasing 
V g brings /oi into resonance with the resonator, reveal¬ 
ing multiple avoided crossings [Fig. 2(b)]. The minimum 
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FIG. 2: Side-gate tuning of a NW circuit and mesoscopic 
Josephson coupling fluctuations, (a) Left inset; false-colored 
optical image of Device 2. Right inset; SEM micrograph show¬ 
ing the single NW junction and the proximal side gate (or¬ 
ange) for voltage control. Sweeping this voltage induces re¬ 
producible fluctuations in the qubit transition frequency /oi. 
Lower inset; example spectroscopy of the qubit transition, 
showing an inhomogeneously broadened linewidth r y/27r = 
13.2±0.3 MHz. A downward trend in /oi is observed as V g de¬ 
creases. At V g < —15 V, foi fluctuates around the resonator 
fundamental (green line), (b) A zoom-in around V g — — 22 V 
shows multiple avoided crossings, (c) At V g — —22.3 V, the 
NW circuit fully hybridizes with the resonator. From the min¬ 
imum splitting, we extract the NW circuit-resonator coupling 
strength g/27v = 34 ± 1 MHz. 


splitting indicates g/2i\ = 34 ± 1 MHz [Fig. 2(c)]. We 
note that while we only perform quasistatic field-effect 
tuning of /oi throughout this experiment, nanosecond 
control should be possible by increasing the bandwidth 
of off-chip filtering. 

We now discuss the impact of charge fluctuations on 
the observed linewidth of the /oi transition, which is of 
interest for qubit applications. Transmon qubits are by 
design insensitive to charge offset fluctuations on the su¬ 
perconducting islands [14], owing to the exponential sup¬ 
pression of charge dispersion when Ej Eq- Field-effect 
control of the Josephson coupling can make /oi sensi¬ 
tive to nearby fluctuating charges. One may expect a 
region with dfoi/dV g = 0 to constitute a charge sweet- 
spot [22] and thus to correlate with a linewidth reduc¬ 
tion. However, we do not observe a correlation between 
the linewidth j/2tt > 10 MHz and \dfoi/dV g \, suggesting 
a different dominant decoherence channel [19]. We sur¬ 
mise a connection between this broad linewidth and the 
soft gap induced in NWs [23] contacted by NbTiN using 
similar fabrication techniques [6]. Parallel experiments 
by the Copenhagen group achieve hard induced gaps in 
epitaxial Al-InAs NWs [24] and ~ 1 fis coherence times 
in a single-JE hybrid transmon [25]. Establishing a clear 
connection between soft induced gap and dissipation in 
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FIG. 3: Flux-bias spectroscopy of a split-junction NW circuit, 
(a) Inset: false-colored optical image showing Device 3 (red), 
its resonator (green), and flux-bias line (purple). Bottom in¬ 
set: SEM micrograph of the two JEs made from one NW. 
Flux-bias spectroscopy shows the tuning of /oi with <f>. (b) A 
high-resolution sweep around <f> = $o/2 shows a strong flux 
dependence of the NW circuit transitions, (c) Measurement of 
resonator transmission around f c with same horizontal range 
as in (b). The avoided crossing of the lowest transition with 
the resonator reveals a much reduced coupling strength. 


Next, we consider a split-junction device where the two 
parallel JEs (each with L = 150 nm) are created from 
one 5 /im long NW [Fig. 3(a)]. As in conventional trans- 
mons, /oi first decreases as flux <f> is threaded through 
the loop. However, near ^ <l>o/2 (4>o = h/2e is the 
flux quantum), a clear departure from transmon-like be¬ 
havior is observed [Fig. 3(b)]. Multiple strongly flux- 
dependent transitions and a new, strong avoided cross¬ 
ing appear symmetrically about 4>o/2. In addition, the 
avoided crossing between the lowest transition and the 
resonator is strongly reduced [Fig. 3(c)] compared to that 
of Device 2 [26]. The highly anharmonic spectrum and 
much reduced coupling to the resonator shows that the 
states associated with the lowest-energy transition are no 
longer plasma modes. 

The observed deviation from the conventional trans¬ 
ition energy spectrum provides a signature of non- 
sinusoidal CPRs in the NW junctions [27]. We now 
show that the observed spectrum can be fully explained 
by the Hamiltonian of a Cooper-pair box (CPB): H = 
4:EqN + V\ (6) + U 2 (2tt<F/$o — 5), provided its split 
junctions do not follow a cosine-shaped Josephson po¬ 
tential [Fig. 4(a)]. Here, the operators N and 5 rep- 



FIG. 4: Fitting theory to the spectrum of the split-junction 
NW circuit, (a) Top inset; equivalent circuit of the device, 
a split Cooper-pair box containing a loop interrupted by two 
NW JEs (phase differences (^ 1 , 2 ), threaded by an externally 
applied flux <f>, and shunted by a capacitance, giving total 
charging energy Ec . Around <f> = 0, the Josephson poten¬ 
tial has a single minimum, producing a weakly anharmonic 
spectrum, (b) The non-sinusoidal current-phase relation of 
the NW JEs determines the particular flux-dependence of the 
transition frequencies around <f> = 4>o/2 as the Josephson po¬ 
tential develops a symmetric double-well profile that tilts as 
<f> is tuned away from 4>o/2. All curves are the result of a 
least-squares non-linear fit of the theoretical model described 
in the text. We identify four fundamental transitions from 
the ground state and three transitions from the first-excited 
state, (c) The CPR calculated from the fit of the spectrum 
to the theoretical model described in the text. 


resent the charge imbalance between islands and the 
phase difference across NW junction 1, respectively. The 
Josephson potential Vi((pi) of junction i is linked to 
its CPR by = (27r/4> 0 ) dVi/dipi, where cp 1 = 5 

and = 27t4>/<Fo — 5. Crucially, we require Vi to be 
27r-periodic but not necessarily cosine shaped. Using 
a simple phenomenological model [19, 21] of the form 

Vi(cpi) = —Kiyjl — Ti sin 2 (ipi /2) and performing a non¬ 
linear least-squares fit with five free parameters, we ob¬ 
tain a quantitative match to all spectral data (best-fit 
values are E c /h = 279 =b 1 MHz, K x /h = 376 ± 13 GHz, 
K 2 /h = 233 ± 2 GHz, T x = 0.86 =b 0.02, and T 2 = 
0.885 =b 0.004). This simple model stricktly applies only 
to short weak links having a single conductance chan¬ 
nel, we therefore use the fit only to extract information 
about the CPRs of the junctions. As shown in Fig. 4(a), 
the corresponding CPRs are evidently skewed. A three- 
parameter fit using Vi corresponding to the CPR of a 
short, diffusive point contact in the many-channel limit 
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[27, 28] showed only slightly worse agreement, as did a 
truncated Fourier series expansion of V{. All approaches 
produce similar skewed CPRs [19]. 

Interestingly, this device can be operated in two dis¬ 
tinct regimes by tuning <f>. Near <F = 0, it operates like 
a transmon, whose eigenstates are plasma modes with a 
weakly anharmonic spectrum. Around <F ~ 4>o/2, it op¬ 
erates like a flux qubit [15] whose two lowest energy levels 
carry opposite persistent currents I p ^ = dEi/d <£>, which 
we estimate to be of order ±100 nA [19]. The possibil¬ 
ity to drive transitions between these distinct persistent- 
current states using coherent microwaves constitutes a 
manifestation of macroscopic quantum coherence [29] in 
our NW circuits. 

In conclusion, we have realized hybrid microwave cir¬ 
cuits made from super-semi NW JEs and characterized 
them using spectroscopy. NW circuits offer several ad¬ 
vantages over traditional aluminum circuits. First, tun¬ 
ing qubit transitions using the field effect in single NW 
JE devices offers an attractive alternative to flux biasing 
of split-junction A1 devices. Second, these NW circuits 
are made exclusively from magnetic-field compatible ma¬ 
terials. Magnetic-field compatible super-semi NW mi¬ 
crowave circuits have the potential to open new avenues 
of research. In particular, very pure solid-state electron 
spin ensembles (e.g., nitrogen impurities in diamond or 
phosphorous donors in silicon) could be field-polarized 
to make coherent quantum memories [30, 31] for hybrid 
quantum processors. In addition, the microwave circuits 
realized here may be useful to control and readout Majo- 
rana bound states [32, 33] in proposed demonstrations of 
non-Abelian exchange statistics [34-37]. Immediate next 
experiments will therefore focus on the study of these 
circuits in up to ~ 0.5 T in-plane magnetic fields. 
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MATERIALS AND METHODS 

The chip was fabricated on a sapphire substrate (single-side polished, C-plane cut, 430 /rm thickness). After cleaning 
the substrate in buffered HF, a NbTiN film (80 nm thickness) was sputtered. Ground planes, resonators and IDCs 
were defined by negative tone electron-beam lithography and reactive ion etching. 

In the next step, In As NWs of typical length 5 — 10 fim and radius 50 — 100 nm were controllably deposited using a 
micromanipulator setup equipped with an optical microscope. Detailed structural and DC transport characterization 
of the InAs NWs used in this experiment have been published in earlier reports [1]. 



FIG. SI: Overview of the chip. 

The fine patterns were designed based on optical alignment to the NWs. The semiconductor-superconductor 
interface area was maximized by covering the NW as much as possible. The overlap between the coarse and the fine 
superconducting structures in each IDC was designed to be several tens of fim 2 . After defining the fine pattern by 
positive tone electron-beam lithography, the semiconductor surface was cleaned in buffered HF for 20 s and then a 
100 nm thick NbTiN film was sputtered in order to overlay the coarse structures as well. In total, twelve devices 
were fabricated on the chip (Fig. SI), and each was coupled to a coplanar waveguide quarter-wave resonator with a 
distinct fundamental frequency in the range 4 to 11 GHz. Table 1 gives an overview of all the devices on the chip. 
Gate voltages in the range ±40 V were applied through on-chip 50 Ut transmission lines with open ends proximal to 
the NW. Off-chip, the lines are filtered by a second-order RC filter and a copper-powder (CuP) filter mounted to the 
mixing chamber plate of the refrigerator (Fig. S2). The two filters have ~ 50 kHz and ^ 1 GHz cutoff frequencies, 
respectively. The JEs of split-junction devices were embedded in a superconducting loop of area 13 x 24 /im 2 . Flux 
through the loop is controlled using on-chip 50 Q transmission lines, with short-circuited termination proximal to the 
loop. Off-chip, the lines are low-pass filtered with nominal 1 GHz cutoff frequency. 
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Table SI. 

Overview of devices on the chip. _Rrt is the two-terminal resistance probed close to the junction at room temperature. 


Device 

Rrt (kfi) 

^design (nm) 

fc (GHz) 

foi (GHz) 

Description 

1 

5.7 

500 

4 

6.7 

single (Fig. 1) 

2 

6.2 

500 

4.5 

6.2 (V e = 0 V) 

single, with gate (Fig. 2) 

3 

2.2 

200 

6.5 

17 (max) 

split, with flux control (Figs. 3&4) 

a 

3.03 

200 

7 

not < 20 

single 

b 

3.07 

200 

6 

13.6 

single, with gate 

c 

short 

50 

11 

- 

single 

d 

open 

100 

8.5 

- 

single, with gate 

e 

short 

50 

10 

- 

single 

f 

3.5 

500 

5 

9.1 (max) 

split, with flux control 

g 

short 

50 

10.5 

- 

single 

h 

1.4 

100 

8 

not < 20 

split, with flux control 

i 

2.6 

100 

9 

not < 20 

single 
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FIG. S2: Wiring schematic of the experimental setup outside and inside the 3 He- 4 He dilution refrigerator (Leiden Cryogenics 
CF-650). Readout (green) and spectroscopic (red) microwave drives are combined in a single coax line (blue) at room tem¬ 
perature. The line is attenuated at the 3 K and 20 mK stages before connecting to the chip feedline. The feedline output is 
isolated from the higher temperature stages by two circulators. The signal is amplified at the 3 K stage by a HEMT amplifier 
(Caltech Cryol-12, 0.06 dB noise figure) and at room temperature by two Miteq amplifiers. The signal is downconverted to 
1 MHz, further amplified, digitized and saved for processing. The setup contains wiring for gate (orange) and flux (purple) 
control on a subset of devices. Both types of control line are low-pass filtered. 
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ADDITIONAL DATA FOR DEVICE 2 



FIG. S3: Reproducibility of the foi fluctuations. Two sweeps of two-tone spectroscopy show the reproducibility of the foi 
fluctuations. These sweeps were separated by two days. 



FIG. S4: Absence of correlation between /oi field-effect sensitivity and linewidth 7 . (A) Qubit transition frequency /01 as a 
function of V g , extracted from Fig. 2. (B) Computed field-effect sensitivity \dfoi/dV g \. (C) Extracted 7 . Vanishing \dfoi/dV g \ 
does not correlate with a reduction in 7 . 
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DATA EXTRACTION FROM FLUX-BIAS SPECTROSCOPY (FIG. 3) 

In order to perform fits to the observed spectroscopy lines, we have extracted a set of data points with error bars 
from the raw data, consisting of a 2D scan of feedline transmission measured as a function of frequency / and voltage 
Vbias applied to the 50 D bias line. We used the following procedure to extract the points [2]. 

1. We applied a Gaussian filter to suppress noise fluctuations, and subtracted the residual background signal. 

2. We converted Wias into flax through the loop, assuming a linear relation 4> = AVb[ as + B. 

3. We identified isolated features in the raw data by removing all points below a transmission threshold. Every 
isolated feature consisted of a connected set of data points surrounding peaks of high transmission in the /-<£> 
plane. 

4. Within every feature, we extract a single data point for each voltage value. The frequency / and uncertainty 
A/ of every point was computed by taking the average frequency of all points in the same feature with the same 
V, weighted by their transmission amplitude. This procedure provided a collection of data points (<£,/, A/). 
The points were manually divided into groups forming continuous /(<£>) transition lines. 


THEORETICAL MODEL AND FITS FOR DEVICE 3 

Many features of the observed transitions coincide with those expected for a double-well potential. The first is the 
linear vanishing of foi at = 4>o/2, which is consistent with the small energy difference between the ground and first- 
excited states in a symmetric double-well potential. The second is the appearance of a strong avoided crossing between 
the two lowest transitions at ~ 0.494> 0 . The avoided crossing naturally arises in a tilted double-well potential, when 
the lowest energy state in the shallower well becomes resonant with the first-excited state of the deeper well. The 
third is the visibility of transitions whose frequency decreases away from <f> = 4>o/2, which is consistent with residual 
thermal population of the first-excited state close to 4>o/2. 

The transition frequencies obtained from the energy levels of the split-junction CPB Hamiltonian, 

H = 4E C N + Vi(<5) + V 2 (27r<f>/$ 0 - 6) . (SI) 

were fit to the extracted data points. The energy levels were computed numerically in the eigenbasis of the charge 
operator TV, truncating the Hilbert space to ~ 200 states such that convergence was reached. Via a Fourier series 
expansion of the 27r-periodic functions 


Vi(lfii) = - X A i,n cos ( n <Pi) , (S2) 

n 

the Josephson terms in the Hamiltonian can be easily expressed in the same basis using the raising and lowering 
operators N± = exp (±T<5), 


Vi(6) = -53^i, n (JV++h.c.), 

n 

F 2 (27T$/$o -S) = -Y,A 2 , n (N + e-^/*°+h.c.). 

n 

Hence, the Fourier series of Vi contains higher harmonics with period 2tt /n, with n > 1, allowing the total Josephson 
energy V\ + V 2 to develop several local minima in a 2ir interval, thereby forming a double-well potential. This situation 
is particularly relevant at <f> ~ 4>o/2, where the odd harmonics of V± and V 2 subtract. Time-reversal symmetry makes 
the minima degenerate at = 4 >q/ 2 (Fig. 4). 

Lacking exact knowledge of the CPR of the junctions, we have tried different phenomenological models for the 

Vito): 


1. The model presented in the text, 


Vito) = -K iy / l-T iS in 2 (^/2), 


(S3) 
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is strictly valid only for short junctions, with total amplitude K{ ~ n^Ao, with Aq the induced gap in the NW, 
resulting from rii transport channels with a typical transparency T^. This model leads to a Hamiltonian with five 
free parameters: Eq , i^i, K 2l Ti, and T' 2 . For the fit we took the Fourier series expansion of V* and truncated 
to seven harmonics in the numerical diagonalization of the Hamiltonian. 

2. The Kulik-Omelyanchuk (K-O) model for the zero-temperature current-phase relation of a diffusive and short 
point contact [3, 4] gives 


nir/2 1 - 

Vi(ipi) = -2Ki J y 1 — sin 2 (^/2) sin 2 xdx. (S4) 

The above relation can be obtained from Eq. (3) by integrating over the probability distribution of transmission 
coefficients T in a diffusive point contact. This model leads to a Hamiltonian with three free parameters, Eq , 
ifi, and K 2 . As before, the Fourier series was truncated to seven harmonics. The best-fit results are shown in 
Fig. S5. 

3. We have also performed a direct fit of the Fourier coefficients increasing the number of coefficients until 

a good agreement with the data was observed. The fit is more sensitive to the Fourier coefficients of the weak 
junction, as the asymmetry in junction strength makes the phase difference across the strong junction stay close 
to zero at any <F. The best-fit values of the seven-parameter fit are Ec/h = 275±1 MHz, Aqi/h = 95.5±0.4 GHz, 
A 2 ,i/ h = 68.9 ±0.2 GHz, A 2A /h = -11.78 ±0.06 GHz, A 2A /h = 3.28 ±0.02 GHz, A 2A /h = -0.98 ± 0.03 GHz, 
and A 2A /h = 0.21 ± 0.03 GHz. Including more Fourier coefficients did not improve the agreement with the 
data, and increased the variance of the best-fit values. 



0/O o 0/O o 

FIG. S5: Fit to the K-0 model. Best-fit values are E c /h = 266 ± 2 MHz, K ± /h = 218 ± 1 GHz, and K 2 /h = 141 ± 1 GHz. 

Due to the non-linearity of the fitting model and correlated errors in the extracted data points, it is not straight¬ 
forward to extract a figure of merit for quantitative comparison of the models. All models require at least one 
non-sinusoidal CPR to fully explain all the features of the data set. 

Although the data are well fit by non-sinusoidal CPRs, there are other mechanisms that can produce similar spectra 
in superconducting circuits. First, the presence of a superinductance is known to produce a potential landscape with 
several minima, as for the fluxonium [5]. NbTiN is known for its high kinetic inductance ±k- Indeed, in our circuit 
we estimate ±k ~ 0.2 nH and loop geometric inductance L geo & 30 pH. The associated total inductive energy 
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FIG. S6: CPRs for the three models. (A) Phenomenological model used in the main text. (B) K-0 model , (C) Fourier 
series using one harmonic for JE 1 and five harmonics for JE 2. 


El = (<3>o/27r) 2 /Z/k ~ hx 800 GHz is, however, much larger than Ej « (K\T\ +K2^2)/4 = h x 132 GHz and therefore 
cannot produce multiple minima (fluxonium requires El < Ej). Second, microwave-induced excitation of Andreev 
bound states (ABS) confined to the junction can also produce a qualitatively similar spectrum [6]. However, ABS 
transitions below 2 GHz require almost perfectly transparent transport channels, which has proven challenging in 
even the most ideal atomic point contacts [6]. We therefore surmise that the JEs remain in their ABS ground state 
and that the observed spectrum is fully due to plasma modes. The potential to observe ABS transitions provides an 
exciting prospect for future work. 


ESTIMATION OF DEVICE 3 PARAMETERS FROM MODEL 

It is possible to estimate interesting properties of Device 3 using the best-fit model values presented in the main text. 
In particular, we can compute the dipolar couplings of the NW circuit to its resonator. The coupling strength of the 
dipole-induced transition between two eigenstates | k) and \l) of the split-junction CPB Hamiltonian is proportional 
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to the corresponding matrix element of TV, n^i = (k\N \l). 

Fig. S7 shows the calculated relative couplings |rzoi|, |no 2 1, |rio 3 1 ? and |rii 2 1 near <F = <Fo/2. Their flux dependence 
explains some interesting features of the data, in particular: 

1. The vanishing visibility of the 0—^2 transition at <£> = 3>o/2. 

2. The strong visibility of the 0^3 transition near <F = <Fo/2. 

3. The reduction of |noi| near Tq/ 2 (see Fig. 3) compared to that of Device 2. 



(0/0 o ) 


FIG. S7: Flux-bias dependence of the NW-circuit-resonator coupling matrix elements for Device 3. The absolute matrix 
element |(/c| TV |Z)| determines the visibility of the transitions in Fig. 3 of the main text. As expected for transmons, only the 

matrix elements \(k\ TV \ k + 1) remain strong away from <f> = <f>o/2. 
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( 0/00 ) 


FIG. S8: Energy levels of the split-junction NW circuit around <f>o/2. (A) Calculated energy levels for the parameters of 
Fig. 4 of the main text. (B) The states belonging to the lowest energy levels, Eq and E \, carry opposite persistent currents 
/ Pj i s= dEi/d&. The two states hybridize at <f> = <f>o/2. Driving /oi with microwaves induces transitions between these two 
macroscopically distinct current states. 
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